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Abstract: We obtain the band edge eigenstates and the mid-band states for the complex, PT-invariant 
generalized associated Lame potentials V PT (x) = —a(a + l)msn 2 (y,m) — b(b + l)msn 2 (y + K(m),m) — 
/(/ + l)msn 2 (y + K(m) + iK'(m),m) — <?(<? + l)msn 2 (y + iK'(m),m), where y = ix+(3, and there are four 
parameters a, b, f,g. This work is a substantial generalization of previous work with the associated Lame 
potentials V(x) = a(a + l)msn 2 (x, m) + b(b + l)msn 2 (x + K(m),m) and their corresponding PT-invariant 
counterparts V PT (x) = —V(ix + (3), both of which involving just two parameters a, b. We show that for 
many integer values of a,b,f,g, the PT-invariant potentials V PT (x) are periodic problems with a finite 
number of band gaps. Further, using supersymmetry, we construct several additional, new, complex, PT- 
invariant, periodic potentials with a finite number of band gaps. We also point out the intimate connection 
between the above generalized associated Lame potential problem and Heun's differential equation. 
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1 Introduction 



In the past few years, Bender and others ^12] have looked at several complex potentials with PT-symmetry 
and have shown that the energy eigenvalues are real when PT-symmetry is unbroken, whereas they occur in 
complex conjugate pairs when PT-symmetry is spontaneously broken. However, there have been relatively 
few papers discussing periodic potentials with PT-symmetry [3J 0]. Recently, we |3J have constructed 
several new classes of analytically solvable, complex, PT-invariant, periodic potentials with the special 
feature that they possess just a finite number of band gaps. The purpose of this paper is to substantially 
increase this list of solvable periodic potentials. 

A few years ago, we obtained the band edges of the associated Lame (AL) potentials [5| 



Here, sn(x,m), cn(x, m), dn(x,m) are Jacobi elliptic functions with elliptic modulus parameter m (0 < 
m < 1). They are doubly periodic functions with periods [AK (m), i2K'(rri)], [4K(m), 2K (m) + i2K'(m)], 
[2K(m),iAK' (m)] respectively |13) . where K(m) = Jq d0[l — m sin 2 #]~ 1//2 denotes the complete elliptic 
integral of the first kind, and K'(m) = K{\ — m). For simplicity, from now on, we will not explicitly 
display the modulus parameter m as an argument of Jacobi elliptic functions. It was shown that the AL 
potentials with integral values of a, b are periodic potentials with a finite number of band gaps We 
also constructed and studied the PT-invariant potentials V PT (x) = —V(ix + (3) obtained from the AL 
potentials via the anti-isospectral transformation of variables x — > ix + (3 [IJ. 

In this paper, we make a substantial generalization of our previous work. We consider the four parameter 
family of generalized associated Lame (GAL) potentials 

V(x) = a(a + l)msn 2 (x, m) + b(b + l)msn 2 (x + K(m),m) 



V{x) 




(1) 



+ /(/ + l)msn 2 (x + K{m) + iK' (m),m) + g(g + l)msn 2 (x + iK' (m),m) 



a(a + l)msn 2 (x, m) + b(b + l)m 2 ' + /(/ + 1) 2 ' + g(g + 1) 

dn (x, m) cn z {x,m) 




2 (x, m) 
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In contrast to the AL potentials of eq. Q where there are two parameters a, 6 and the two terms correspond 
to real translations of the independent variable x by and K(m), the GAL potentials of eq. @ have four 
parameters a,b,f,g and the four terms correspond to complex translations of the independent variable x 
by 0, K(m), K(m) +iK'(m),iK'(m). Although the GAL potentials are real, they do have singularities on 
the real axis coming from the zeros of the Jacobi elliptic functions sn (x) and cn (x) in the last two terms. 
Consequently, we will focus our attention on the PT-invariant versions of the GAL potentials, which are 
given by 

V PT (x) = —a(a + l)msn 2 (y, m) — b(b + l)msn 2 (y + K(m),m) 

— /(/ + l)msn 2 (y + K{m) + iK'(m),m) — g(g + l)msn 2 (y + iK' '(m),m) 

I i i\ 2/ n ? /i I i \ cn 2 (y,m) dn 2 (y,m) . 1 

= -a{a + ljmsn (y, m) - bib + l)m— - /(/ + 1) — T . - g{g + 1)- 



where 



dn 2 (y,m) cn 2 (y,m) sn 2 (y, m) 

[a(a + l), 6(6 + !),/(/ + !), g(g + l)], (3) 



y = ix + (3, (4) 



with j3 being an arbitrary constant. We shall frequently use the notation [a(a + l), 6(6+1), f(f + 1), g(g + l)] 
to denote V PT {x). In this notation, PT-invariant ordinary Lame potentials are denoted by [a(a+l), 0, 0, 0], 
and PT-invariant AL potentials are denoted by [a(a + 1), 6(6 + 1),0,0]. Here, the arbitrary constant j3 
is chosen so as to avoid the singularities of the Jacobi elliptic functions on the real axis. We show that 
several of these periodic potentials for specific integer values of a, 6, /, g have a finite number of band gaps. 
Looking at the symmetry of these potentials, we are in fact tempted to conjecture that many (and perhaps 
all) of these potentials with integral values of the parameters a, 6, /, g also have a finite number of band 
gaps. It would be nice if this conjecture could be proved. 

In addition, we also discover a huge class of mid-band states when at least one of the parameters a, 6, /, g 
is half integral. As a special case, we find some new mid-band eigenstates of the associated Lame potentials. 
Further, we show that the Schrodinger equation for the generalized AL potential is intimately connected 
with the celebrated Heun's differential equation [7j. In fact, using the exact solutions obtained in this paper, 
one can immediately obtain the corresponding solutions of Heun's equation. In another related paper jS], 



we use this connection and discover a wide class of new quasi-periodic solutions of Heun's equation. 

Finally, using the exact eigenstates of the GAL potentials Q and the machinery of super symmetric 
quantum mechanics [S], we construct several more potentials with finite band-gaps. There is one important 
point involved here using which we are in fact able to construct many more supersymmetric partner 
potentials corresponding to a given potential. The key point to note is that normally, in supersymmetric 
quantum mechanics 0, given a potential V-(x), the ground state wave function i/jq(x) is used to construct 
the superpotential W^x) = —^' (x)/ipo(x), which then yields the supersymmetric (SUSY) partner potential 
V+(x) = W 2 + W . If one uses any excited state wave function ifi(x) of V-(x) to construct a superpotential 
W(x), then the original potential V-{x) is recovered correctly (by construction), but the corresponding 
partner potential U+(x) turns out to be singular on the real x-axis due to the zeros of the excited state 
wave function ip{x). However, as has been noticed recently |1U| . if we consider PT-symmetric complex 
potentials, then the singularity is not on the real axis. Besides, as we have stressed previously [H 111), 
in the case of doubly periodic potentials composed of Jacobi elliptic functions, both V{x) and V PT {x) 
can be simultaneously periodic even though their periods are different. In this way, by starting from the 
analytically solvable Lame and associated Lame potentials and using the excited state band edges of the 
corresponding PT-symmetric potentials, we discover a wide range of new, analytically solvable, complex 
PT-invariant periodic potentials with a finite number of band gaps. As an illustration, we discuss a few of 
these potentials in detail. 

The plan of the paper is the following. In Sec. 2 we discuss the PT-invariant GAL potentials (J3J) in some 
detail and obtain band edges as well as mid-band states of several of these potentials. As a byproduct we 
also obtain some new solutions of the AL potentials (which we had missed in earlier work [BJ ) . Further, 
we show that the class of potentials [a{a + 1), 0, 0, g{g + 1)] have finite number of band-gaps in case a, g are 
integers. In Sec. 3 we start from the energy eigenstates obtained in Sec. 2 and using both the ground state 
as well as excited state wave functions, obtain new periodic PT-invariant potentials with a finite number 
of band gaps. In Sec. 4 we briefly discuss the connection between the solutions of the potentials © and 
Heun's differential equation. 
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2 Solutions for the Generalized Associated Lame (GAL) Potentials 

A few years ago, we obtained analytic solutions of the associated Lame potentials Q 5,6 and showed 
that when a, b are integers, then the resulting potentials all had a finite number of band gaps. The purpose 
of this section is to show that the complex PT-invariant GAL potentials as given by eq. © are also quasi- 
exactly solvable. In particular, we show that the band edges or mid-band states of these problems can be 
obtained depending on whether a + b + / + g (or a — b — f — g) is an integer or an arbitrary non-integer 
number. It should be noted that we are considering PT-invariant potentials ©, since the corresponding 
real potentials (J2J) are singular on the real axis. 

It may be worthwhile explaining the underlying basic idea here, even though it has been well established 
by us before [I]. Note that if ip{x) is a solution of the Schrodinger equation for the real potential V{x) with 
energy E, then tp^ix + (3) is a solution of the Schrodinger equation for the complex potential —V(ix + (3) 
with energy — E, where (3 is an arbitrary nonzero constant. The new potential — V(ix + /?), generated 
by the anti-isospectral transformation x — > i x + (3 |12) . is clearly PT-symmetric and will be denoted by 
V PT {x). Further, if ip(x) and ip(ix + (3) satisfy appropriate boundary conditions, they are eigenfunctions 
of V(x) and V PT (x) respectively. The ordering of energy levels for V PT {x) is the opposite of the ordering 
of energy levels for V(x). 

In this paper, our main focus is on the Schrodinger equation (h = 2m = 1) 

d 2 

~ tt^(^) + V PT {x)^{x) = Eijj(x) , (5) 
ax z 

where V PT (x) is the potential given by eq. Eq. © is called the generalized associated Lame equation, 
and we are seeking its eigenstates and mid-band states. 

2.1 Symmetries 

At this stage, it is worth pointing out the symmetries of the PT-invariant GAL potential © and hence 
the corresponding Schrodinger equation (jSJ. 

1. The potential (© and hence the Schrodinger eq. © remains unchanged when any one (or more) of 
the four parameters a, b, /, g change to —a — 1, —b — 1, — / — 1, —g — 1 respectively. 
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2. Under the translation y — > y + K(m), the GAL potential [a(a + 1), 6(6 + 1), /(/ + 1), y(y + 1)] goes to 
the potential [b(b + 1), a(a + l),g(g + 1), /(/ + 1)]. Hence, both GAL potentials must have the same 
energy eigenvalues and the corresponding energy eigenfunctions are simply related: y — > y + K(m), 
i.e. 

E PT (b, a, g, /; m) = E PT (a, b, /, y; to) , ip(y, b, a, g, /; m) oc ip(y + K(m),a, 6, /, g; m) . (6) 

3. Similarly, by considering the translations y — > y + K{m) + iK'(m), and y — > y + iK'(rn), it is easy 



E PT (f, g, a, b; m) = E PT (a, b, f, g; m) , ^(y, /, y, a, 6; m) oc V>(y + K(m) + iK'(m),a, b, /, y; m) . (7) 

E PT (g, f, b, a; m) = E PT (a, b, f, y; m) , ip(y, g, f, b, a; m) oc i/)(y + iK'{m),a, b, f, y; m) . (8) 

Thus, once we obtain the eigenvalues and eigenfunctions of a given GAL potential [a(a + 1), b(b + 1), /(/ + 
l))^!^ + !)]> then we immediately know the eigenvalues and eigenfunctions of three other potentials: 



[b(b+l), o(o+l), y(y+l), /(/+!)], [/(/+!), y(y+l), a(o+l), 6(6+1)] and [y(y+l), 6(6+1), a(o+l)]. 



Therefore, it suffices to present results for only one of the four potentials. 
2.2 Duality Relations 

We shall now derive some remarkable relations relating the quasi-exactly solvable eigenvalues and eigen- 
functions (corresponding either to the band edges or mid-band states) of two GAL potentials at two 
different values m and 1 — m of the modulus parameter. 

To that purpose we start from the Schrodinger eq. (JSJ) for the PT-invariant GAL potential ((31) • On 
using the relations ^3 



to show that 



\J~m sn(y, m) 



dn [iy + K (m) + iK(m), 1 — m] 



dn(y, m) = \/T 



msn [iy + K'(m) + iK(rn), 1 — m] 



^/mcn(y,rn) = iVT 



men [iy + K'(m) + iK(m), 1 — m] , 



(9) 
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and denning a new variable w = iy + K'(m) + iK(m), the Schrodinger eq. (JSJ) takes the form 

r(w) - a a + 11- m sn 2 («;, 1 - m + 5 g + 11 - m)-^- >-f(f + 1 — aHi ( 

dn (u>, 1 — m) cn z (u>,l — mj 



+6(6 + 1) — 27 — ^ Mw) = -[a(a + 1) + 6(6 + 1) + /(/ + 1) + g(g + 1) + E)i>(w) . (10) 

sn z (u>, 1 — m) 

On comparing eqs. © and IjlUjl we then have the remarkable relations 

E PT (a, b, f, g, m) = -[a(a + 1) + 6(6 + 1) + /(/ + 1) + g(g + 1)] - E PT (a, g, f,b,l-m), 

ip(y, m) oc ip(iy + K (m) + iK(m), 1 — m) , (11) 

which is valid for the QES states corresponding to either the band edges or mid-band states. Note that 
here, a, 6, /, g can be arbitrary (real) numbers and are not restricted to integer values. This is a very 
powerful relation which has several interesting consequences. One immediate important consequence of 
eq. Hll|) is that for arbitrary integer values of a, g, the potential [a(a + 1), 0,0, g(g + 1)] has only a finite 
number of band-gaps. This happens because, for f = g = 0, one has 

E PT (a, b, 0, 0, m) = -[a(a + 1) + 6(6 + 1)] - E PT (a, 0, 0, g = b, 1 - m) , (12) 

so that both the potentials must have the same number of band-edges and band-gaps and we have already 
proved [S] that the AL potentials have finite number of band gaps in case a, b are integers. 

2.3 QES Solutions 

Let us now seek solutions of the Schrodinger eq. (JSJ for the PT- invariant GAL potential ©• On making 
the ansatz 

ip(x) = dn- b (y)sn~ 9 (y)cn- f (y)(j)(y) , y = ix + (3 , (13) 

it is easily shown that 4> satisfies the equation 

±ii l s , or ,sn(y)cn(y) cn(y)dn(y) dn(y)sn(y) 2 

(j) (y) + 2[mb — -— g — + / — ]<j> (y) + [Qmsn [y] - R]<f>{y) = , (14) 

dn(y) sn(y) cn(y) 

where 

Q = (b + g + f)(b + g + f-l)-a(a + l), R = E + (/ + g) 2 + m(g + 6) 2 . (15) 

It is well known |14j that this is a quasi-exactly solvable (QES) problem. We shall now systematically 
consider solutions of eq. (|14|) for several special cases and then finally consider the most general case. 



2.4 b = f = g = 

The simplest possibility is when three out of the four parameters a, b, /, g are zero. For example, when 
b = f = g = 0, then the problem reduces to the PT-invariant version of the well studied Lame potential 
problem. We might add here that, instead of a, if any one of the other parameters b, f, g is nonzero, one 
still has a potential which is strictly isospectral to the PT-invariant Lame potential. It may be noted 
that while the Lame potential is a periodic potential with (real) period 2K(m), the PT-invariant Lame 
potential has real period 2K'(m). Further, the band edge eigenvalues, eigenfunctions and the discriminant 
A of V PT {x) are related to those of Lame potential by @] 

E PT (m) = -E 2a -j(m) , ip PT (x,m) oc ip^a-j (ix + /3, m) , j = 0,1,2, ... ,2a 
A PT (E,m) = A[JS + o(o + l),l-m]. (16) 

From eq. I|ll|l. it follows that the PT-invariant Lame band-edge eigenvalues and eigenfunctions, for 
integral a satisfy the remarkable relations (j = 0, 1, 2, ...2a) 

E PT (m) = -a(a + 1) - E^^m) , ipj(y, m) oc Tp 2a -j(w + K'(m) + iK(m), 1 - m) . (17) 

We would like to add here that even the mid-band states satisfy (for half-integral a) relations analogous 
to ltT7|): 

Ej(m) = a(a + 1) - E a _ 1/2 -j(m) , tpj(y,m) oc tp a -i/2-j{w + K ' X™) +iK(m), 1 - m) , (18) 

where j = 0,1, 2,..., a — 1/2 . Note the remarkable fact that for any integer a, all bands and band gaps 
exchange their role as one goes from the Lame potential to its PT-invariant version V PT (x) [1]. 

The next simple possibility is when two of the four parameters o, b, /, g are zero. Here there are three 
distinct possibilities which we discuss one by one. 

2.5 f = g = 

In this case the problem reduces to the PT-invariant AL potential which we have already discussed at great 
length [SJIH1- Note that if either a or b is zero (or -1), then this potential reduces to the PT-invariant of the 
Lame potential. As previously shown by us [H], for arbitrary integral values of a and b, AL potentials are 
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exactly solvable problems with finite number of band-gaps for which one can write down the form of all the 
band edge eigenfunctions, as we do below. We note here that when a > b are both integers, then there are 
precisely a bound bands (some of which are unusual in that both the band edges are of the same period), 
same (a) number of band gaps and all the 2a + 1 band edges are analytically known beyond which there is 
a continuum band extending to E = oo. Note that if b > a, then also there are b bound bands and b band 
gaps and the corresponding eigenfunctions are simply obtained from the a > b case by the transformation 
x — > x + K(m) while the a = b case essentially corresponds to the Lame potential [a(a + 1), 0, 0, 0]. Without 
any loss of generality, we shall only consider AL potentials with a > b. 

The form of the 2a + 1 band edge eigenfunctions of the AL potential depends on whether a — b is an 
odd or an even integer. For example, when b = a — 2p — 1 (p = 0, 1, 2, ...), then there are: 

p eigenstates of the form sn(x)cn(x)dn(x)F p _i [sn 2 (x)], 

p+1 eigenstates of the form dn a_2p (x)-F p [sn 2 (x)], 

a — p eigenstates of the form cn(x)dn 2p+1 ~ a (x)F a _p_i [sn 2 (x)], 

a — p eigenstates of the form sn(x)(dn 2p+1_a (x).F a _p_i[sn 2 (x)]. 
On the other hand, when b = a — 2p, there are: 

p eigenstates of the form cn(x)dn a_2p+1 (x)F p _i[sn 2 (x)], 

p eigenstates of the form sn(x)dn a ~ 2p+1 (x)F p _i [sn 2 (x)], 

a — p eigenstates of the form sn(x)cn(x)dn 2p_a (x)F a _ p _i[sn 2 (x)], 

a — p+1 eigenstates of the form dn 2p_a (x)F a „ p [sn 2 (x)]. 
Here F n [sn 2 (x)] denotes a polynomial in sn 2 (x) of order n. 

We would like to re-state here that all the eigenstates of the PT-invariant version of the AL potentials 
are immediately obtained from the known eigenfunctions of the associated Lame problem and the ordering 
of energy levels of these is the opposite of the corresponding AL problem. Hence, this is also an exactly 
solvable problem with a finite number (a) of band gaps and 2a + 1 known band edges when both a, b are 
integers. 
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2.6 b = f = 

Following our discussion for the AL case, without any loss of generality we assume here that a > g. In this 
case, one obtains n + 1 QES solutions when a + g = n (or g — a = n + 1) with n = 0, 1, 2, .... The QES 
solutions for n = 0, 1, 2, 3, 4 are given in Table 1. In particular, for any choice of a(a + 1), Table 1 lists the 
eigenstates for various values of g(g + 1). The general form of these eigenfunctions is obtained from the 
corresponding AL eigenfunctions as given in Table 3 of [5] by simply interchanging dn(y) and sn(y). 
A few remarks are in order. 

1. Since we are considering the case (b = f = 0), the duality relation (jlljl takes the form 

E PT (a, 0, 0, g; m) = -[a(a + 1) + g(g + 1)] - E PT (a, b = g, 0, 0; 1 - m) , 

ip(y, a, 0, 0, g; m) oc ip(iy + K'(m) + iK{m), a, b = g, 0, 0; 1 — m) , (19) 

Using Table 3 of ref. [5] and this duality relation, it is straightforward to obtain all the QES 
eigenstates, thereby providing an independent check on the results given in Table 1. Further, it 
follows that for arbitrary integer values of a and g, [a(a + 1), 0,0, g(g + 1)] is an exactly solvable 
potential problem with a finite number (a) of band-gaps. From the duality relation ()19|) . it follows 
that for integer values of a, g 

E PT (a, 0, 0, g; m) = - [a(a + 1) + g(g + 1)] + E 3 (a, b = g, 0, 0; 1 - m) , (20) 

and hence the corresponding discriminants A are related by 

A PT (E, m; a, 0, 0, g) = A[E + a(a + 1) + + 1), 1 - m; a, b = g, 0, 0] . (21) 

2. Following the structure of the eigenfunctions of the AL potentials as given above, it is now straight- 
forward to write down the general form of the eigenfunctions for arbitrary value of n. However, to 
obtain the corresponding eigenvalues, one needs to solve cubic and higher order equations. 

3. Under the transformation y — > y + iK'{m) followed by the interchange of a and g (note b = f = 0), 
the Schrodinger eq. © for the GAL potential (j3J) remains unchanged. Thus it follows that under 
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the interchange of a with g, the eigenvalue spectrum must remain unaltered. Clearly, this is only 
possible if either the energy eigenvalues remain unchanged under this transformation, or if two of the 
eigenvalues go into each other. It is easy to verify from Table 1 that the eigenvalues corresponding 
to the eigenfunctions of period 2iK'(m) remain unaltered under a —* g while the other eigenvalues 
go into each other under this transformation. 

4. Similarly, From Table 3 of 5 , it is easy to check that for the AL potentials l(T|). the eigenvalues 
corresponding to the eigenfunctions of period 2K(m) remain unaltered under a — > b while the other 
eigenvalues go into each other under this transformation. This happens because the AL potentials 
remain unaltered under the transformation y — ► y + K(m) followed by the interchange of a with b. 

Summarizing, we have discovered new exactly solvable potential problems with a finite number of 
band gaps when a, g are arbitrary integers. In fact everything about these potentials can be derived from 
previously known results for AL potentials. 

2.7 b = g = 

In this case, one obtains n + 1 QES solutions when a + f = n with n = 0, 1,2, .... The solutions for 
n = 0, 1, 2, 3, 4 are given in Table 2. In particular, for any choice of a(a + 1), Table 2 lists the eigenstates 
for various values of /(/ + 1). The general form of these eigenfunctions is simply obtained from the 
corresponding AL eigenfunctions as given in Table 2 of [3] by interchanging dn(y) and cn(y). 
Some comments are in order at this stage. 

1. The form of eigenfunctions for arbitrary value of n is easily written down following the structure of 
the AL eigenfunctions given in the last section. 

2. From eq. (|1 lj) it follows that the potential © with b = g = is a self-dual potential, satisfying 

E? T (a, /, m) = -[a(a + 1) + /(/ + 1)] - ^ T (a, /, 1 - m) . (22) 

Using Table 2, it is easily checked that indeed this is true, for any values of a, /. In particular, 
whereas 5§, <5g are invariant under m — > 1 — m, 5q <-> 5j under the same transformation. 
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3. Under the transformation y — > y + K{m) + iK'{m) followed by the interchange of a and / (note 
b = g = 0), the Schrodinger eq. © for the GAL potential © remains unchanged. Thus it follows 
that under the interchange of a with /, the eigenvalue spectrum must remain unaltered. Clearly, 
this is only possible if either the energy eigenvalues remain unchanged under this transformation, or 
if two of the eigenvalues go into each other. It is easy to verify from Table 2 that the eigenvalues 
corresponding to the eigenfunctions of period 2K(m) + 2iK'(m) remain unaltered under a — > / while 
the other eigenvalues go into each other under this transformation. In particular, while 65,6s are 
invariant under a f , 6q <-> 5j under the same transformation. 

2.8 / = 

Let us consider the case when only one out of the four parameters a, b, /, g is zero. As an illustration, we 
discuss the case / = 0. In fact, as described below, once we know the eigenstates of this problem, the 
eigenstates of the other three problems corresponding to either b or g or a equal to zero are immediately 
obtainable, since the four potentials are related by translations of the independent variable. 

For the case / = 0, one obtains ^^{^2^) QES solutions when n is even (odd). Here a + b + g = n 
with n = 0, 1, 2, .... The QES solutions for n = 0, 1, 2, 3 are given in Table 3. In particular, for any choice 
of a(a + 1), Table 3 lists the eigenstates for various values of (b + g)(b + g + 1). 

Some remarks are appropriate. 

1. By looking at the structure of the QES eigenfunctions in Table 3, it is easy to write down the nature 
of eigenfunctions for the general case. 

2. From Table 3, it is easily checked that the duality relation 

E PT (a, b, g,m) =- [a(a + 1) + b(b + 1) + g{g + 1)] - E PT (a, g,b,l-m). (23) 
is indeed satisfied. In particular, both 69, 5±o are invariant under b <-> g followed by m — > 1 — m. 

3. Under the transformation y — > y + K(m) followed by the interchange of a and b, and replacing g by /, 
the Schrodinger eq. © for the GAL potential (j5J) with / = goes over to the Schrodinger equation 



12 



for the GAL potential © with g = 0. Hence, under the interchange of a and b and replacing g by /, 
all the energy eigenvalues of the potential © with / = must go over into those of © with g = 0, 
while the corresponding eigenfunctions are simply obtained from Table 3 by replacing y by y + K{m). 

4. Using similar reasoning it also follows that under the interchange of a with g and replacing b by /, 
all the energy eigenvalues of the GAL potential © with / = go over to those of potential © with 
6 = while the corresponding eigenfunctions are obtained from Table 3 by replacing y by y + iK'(m). 
And finally, under the interchange of b with g and replacing a by /, all the energy eigenvalues of the 
GAL potential © with / = go over to those of potential © with a = 0, while the corresponding 
eigenfunctions are easily obtained from Table 3 by replacing y by y + K(m) + iK'(m). 

2.9 The General Case: a, b, f, g All Nonzero 

Finally, let us discuss the most general case when all the four parameters are nonzero. In this case one 
obtains n + 1 solutions when a + b + / + g = 2n with n = 0, 1, 2, .... The QES solutions for n = 0, 1 are 
given in Table 4. 

1. It is easy to see that in the general case, the eigenfunction is of the form 

n 

i> = sn-3(y)cn^(y)dn- b (y) ]T A k sn 2k (y) , (24) 

fc=0 

while the corresponding eigenvalues are solutions of a n + l'th order equation. 

2. It can be checked from Table 4 that 5\\ is invariant under b *-* g followed by m — > 1 — m. 

3. The GAL potential © and hence the corresponding Schrodinger eq. © is invariant under the 
transformation y — > y + K(m) followed by the interchange of o with b and / with g. Hence, under 
the interchange of a with b and / with g, all the eigenvalues of the GAL system must either remain 
invariant or go into each other. In fact it is easily checked from Table 4 that all the eigenvalues are 
invariant under the interchange of a with b and / with g. Extending this argument, in fact one finds 
that all the eigenvalues are also invariant under a f,b «-> g as well as under a «-> g,b <-> / . 
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2.10 Mid-Band States 

So far we have discussed the results for the PT-invariant GAL potentials, which give eigenvalues and 
eigenfunctions corresponding to the band edges. It may be noted that in all these cases, while a,b,f,g 
need not be integers, either a + b + f + g or a — b — f — g is always integral. We now show that when at 
least one of a, b, /, g is half- integral and either a + b + f + g and/or a — b — f — y is an arbitrary number 
(being an integer is of course a very special case here), then one can obtain doubly degenerate eigenstates 
which correspond to mid-band states. In fact depending on whether we want b or / or g to be half-integral 
(with the other two parameters being integral), we need to use different trial solutions. Therefore, we shall 
consider all three cases one by one. 
Case 1: b half- integral 

We start from eq. (|14[) and further substitute the ansatz 

<P(y) = [cn(y)+isn(y)] t Z(y), (25) 

where t is any real number. After lengthy but straightforward algebra, one can show that Z(y) satisfies 
the equation 

Z»(y) + [m M y) + 2rnb S -^^l - 2g^^l + 2f^^]Z>(y) 

dn{y) sn(y) cn(y) 



+ [-{R + t z ) + {Q + t z )msn\y) - 2itg^{ + 2it/(l - m)- 

sn(y) cn(y) 

+imt(2b +2f + 2g- l)sn(y)cn(y)]Z(y) = , (26) 

where R and Q are as given by eq. (|15|) . Not surprisingly, Z(y) =constant is a solution with energy 
E = -(4t 2 + m)/4 provided f = g = 0,b= 1/2, a = t - 1/2 (i.e. b + f + g = 1/2). 

One can build solutions for higher values of b + f + g from here. In particular, for b + f + g = 2M + 1/2, 

we consider the ansatz (M = 0, 1, 2, ...) 

M M-l 
Z(y) = £ A k sn 2k (y) + cn(y)sn(y) £ ^sn 2fe (y) , (27) 
fc=0 fc=o 

while if b + f + g = 2M + 3/2 then we consider the ansatz (M = 0, 1, 2, ...) 

M M 
Z(y) = cn(y) £ A k sn 2k (y) + sn(y) £ £ fc sn 2fe (y) , (28) 

A;=0 fc=0 
14 



Substitution into eq. (|26[) and simplification yields analytic expressions for the energy eigenvalues and 
eigenfunctions for arbitrary M for 6 = 1/2 and 6 = 3/2. In particular, for 6 = 1/2, we find that 

b = l/2,f = p,f + g = N,a = t-l/2, E = -[t 2 + m(g + 6) 2 ] , (29) 

where both /, g are nonnegative integers satisfying / + g = N with N = 0,1, 2, .... 
Similarly, when 6 = 3/2, a = t — 1 /2, / = p, f + g = N we find that 

£ = m(2g + 1) - [1 + t 2 + m(s + 6) 2 ] ± ^ {2g + l) 2 m 2 + 4m(iV + 1)(/ - g) + 4(1 - m)t 2 (30) 

where, / and g are again nonnegative integers. In all these cases, the corresponding eigenfunctions have the 
form as given above in eqs. (|27|) and (|28[). For small values of N, the explicit coefficients ^4^, appearing 
in the eigenfunction expressions can be easily written down. For example, for 6=1/2 and N = 1, the 
eigenfunction is Z(y) = Acn(y) + Bsn(y) with = it in case / = 1,5 = while § = i in case g = 1, / = 0. 

For the special case of / = g = and t 7^ 1/2, these results represent the generalization of results 
obtained by us previously [0] in the case of AL potential. Further, for / = g = 0, t = 1/2, the results 
obtained above match with the energy eigenvalue expressions obtained in ref. [H] (as they should!). 

Several comments can be readily made. 

1. Since, in the variable y, the GAL potential © has period 2K(m) as well as 2iK'(m), hence ipiu) 
and ip{y + 2K(m)) as well as ip{y + 2iK'(m)) are all eigenfunctions of GAL equation with the same 
eigenvalue. As a consequence, 4>(y) = [cn(y) — isn(y)]*Z(y) is also the eigenfunction with the same 
eigenvalue. Thus for any nonintegral t, each level is doubly-degenerate. The same remark also applies 
to the other two solutions (when / or g is half integral) discussed below. 

2. There is one remarkable symmetry associated with eq. (|26|). In particular, notice that this equation 
is invariant under t — ► — t followed by i — ► — i (where i = y/— 1). But under this transformation, the 
ansatz becomes 

<p(y) = [cn(y, m) - isn(y, m)p' , (31) 

Hence it follows that the energy eigenvalues must be independent of sign of t, i.e. they must be a 
function of t 2 . Similar remarks also apply in the other two cases discussed below (i.e. when /, g are 
half-integral) . 
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3. For integral t, both a, b are half integral and these solutions reduce to those discussed in the last 
section and in that case they correspond to QES band edge eigenstates. 

4. Here we have obtained solutions ip{y) in which a = t — 1/2, f = p, g = N — p and 6=1/2 or 3/2. 
In view of the symmetries of the GAL potentials, we then also have solutions ip{y + K(m)) with 
the same energy in case b = t — 1/2, g = p, f = N — p and a is either 1/2 or 3/2. Similarly we 
have solutions ip(y + K(m) + iK'(m)) with the same energy in case f = t — 1/2, a = p,b = N — p 
and g = 1/2 or 3/2. Further, we also have solutions ip(y + iK'(m)) with the same energy in case 
g = t- 1/2, a = N-p,b = p and / = 1/2 or 3/2. 

Case 2: / half-integral 

We start from eq. (|14[) and further substitute the ansatz 

<P(y) = [dn(y)+ik S n(y)] t Z(y), (32) 

where t is any real number and k = \fm. After some lengthy but straightforward algebra, one finds that 
Z(y) satisfies the equation 

v" i \ . ro-w r \ i o u sn (y)™(y) cn(y)dn(y) dn(y)sn(y) , 

Z (y) + [2iktcn(y) + 2mb — — 2g + 2/ — \Z (y) 

dn(y) sn(y) cn(y) 

+[-(R + mt 2 ) + (Q + t 2 )msn 2 (y) - 2itkg^\ - 2iktb{\ - m)^\ 

sn(y) dn(y) 

+ikt(2b + 2f + 2g- l)sn(y)dn(y)]Z(y) = , (33) 

where R and Q are as given by eq. (fT3|) . Not surprisingly, Z(y) =constant is a solution with energy 
E = -(4mt 2 + l)/4 provided b = g = 0, / = 1/2, a = t- 1/2 (i.e. 6 + / + 5 = 1/2). 

One can build solutions for higher values of b+f+g from here. In particular, in case b+f+g = 2M+1/2, 
we consider the ansatz (M = 0, 1, 2, ...) 

M M-l 
Z(y) = £ A fc sn 2fc (y) + sn(y)dn(y) £ 5 fc sn 2fc (y) , (34) 

k=0 k=0 

while if b + f + g = 2M + 3/2 then we consider the ansatz (M = 0, 1, 2, ...) 

M M 
Z{y) = dn(y) £ A fc sn 2fc (y) + sn(y) ^ 5 fc sn 2fc (y) , (35) 

fc=0 k=0 
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On substituting this ansatz in eq. ()33|) and making algebraic simplifications, we obtain analytic expressions 
for the energy eigenvalues and eigenfunctions for arbitrary M for / = 1/2 and / = 3/2. In particular, for 
/ = 1/2, we find that 

f = 1/2, b + g = N, a = t-l/2, E = —[mt 2 + (y + f) 2 ] , (36) 

where both b, g are nonnegative integers satisfying b + g = N with N = 0,1, 2, .... 
Similarly, when / = 3/2, a = t — l/2,b + g = N we find that 



E = (2g + 1) - [(1 + t 2 )m + (y + f) 2 } ± yj (2g + l) 2 + Am(N + 1)(/ - g) - 4m(l - m)t 2 (37) 

where, b and g are again nonnegative integers. In all these cases, the corresponding eigenfunctions have 
the form as given above in eqs. (|34|) and (|35|) . For small values of N, the explicit coefficients Ak,B^ in the 
eigenfunction expressions can be easily written down. Further, as in the half-integral b case, one can write 
down three more solutions with the same energy. 
Case 3: g half- integral 

We start from eq. ()14j) and further substitute the ansatz 

<f>(y) = [dn(y) + fccn(y)]'Z(y) , (38) 

where t is any real number. After algebraic simplification, it is easy to show that Z(y) satisfies the equation 

7"( \, r out f\±o u sn (y) cn (y) o cn (y)Mv) . OJ dn(y)sn(y) ; 
z W + l-2ktsn(y) + 2mb 2g — + 2/ — — }Z (y) 

Mv) sn (y) cn (y) 

+[-R +(Q + t 2 )msn 2 (y) - 2ktb^\ - 2kt /^M 

dn(y) cn(y) 

+H(26 + 2/ + 2g - l)cn(y)dn(y)]Z(y) = , (39) 

where i? and Q are as given by eq. (|15[). Not surprisingly, Z(y) = constant is a solution with energy 
E = -(1 + m )/4 provided b = f = 0,g = 1/2, a = t - 1/2 (i.e. b + f + g = 1/2). 

One can build solutions for higher values of b+f+g from here. In particular, in case b+f+g = 2M+1/2, 
we consider the ansatz (M = 0, 1, 2, ...) 

Af A/-1 
Z(y) = ]T A fc sn 2fc (y) + cn(y)dn(y) ^ S fe sn 2fe (y) , (40) 
fc=0 fc=o 
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while if 6 + / + g = 2M + 3/2 then we consider the ansatz (M = 0, 1, 2, ...) 

M M 
Z(y) = cn(y, m) £ ^ fc sn 2fc (y) + dn(y) £ i? fc sn 2fc (y) , (41) 
fe=0 fc=o 

Substituting this ansatz in eq. (|39j) and simplifying, one gets analytic expressions for the energy eigenvalues 
and eigenfunctions for arbitrary M for 6 = 1/2 and b = 3/2. In particular, for b = 1/2, we find that 

<7 = l/2, b + f = N, o = t-l/2, £ = -[(/ + 5 ) 2 + m( 5 + 6) 2 ], (42) 

where both 6, / are nonnegative integers satisfying b + f = N with iV = 0, 1, 2, .... 
Similarly, when 5 = 3/2, a = t — 1/2, b + f = N we find that 

£ = 1 + 2/ + (26 + l)m - [(/ + gf + m(g + b) 2 ] ± yj(l-m) [(2/ + l) 2 - (26 + l) 2 m] + 4mt 2 (43) 

where, 6 and / are again nonnegative integers. In all these cases, the corresponding eigenfunctions have 
the form as given above in eqs. (|4U|) and (|41|). For small values of N, the coefficients ^^,5^ appearing in 
the eigenfunctions can be easily written down. Further, as in the half-integral 6 case, one can write down 
three more solutions with the same energy. 

3 Supersymmetry and Potentials with a Finite Number of Band Gaps 

We shall now start with the ground state as well as the excited state eigenfunctions of various PT-invariant 
GAL potentials discussed in the last section and using supersymmetry obtain the corresponding SUSY 
partner potentials. In this manner, we obtain many new periodic potentials V+ (x) with a finite number of 
band gaps. As emphasized in the introduction, unlike real potentials, if we take a complex PT-invariant 
potential, then even if we start with an excited state wave function and calculate the corresponding 
superpotential W, the singularities in W and hence V+(x) are not on the real axis, and do not cause 
problems. 
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3.1 Supersymmetry Partners of PT-Invariant Lame Potentials 

The simplest case is when only one parameter, (say a) is nonzero. This gives the PT- invariant Lame 
potential 

V(x) = -a(a + l)msn 2 (y) . (44) 

For concreteness, take a = 1, which yields V(x) = — 2msn 2 (y). Here, the three band edge eigenfunctions (in 
order of increasing energy eigenvalues) are sn(y),cn(y), dn(y). It is easily computed that corresponding to 
these three eigenstates, the corresponding partner potentials (up to a constant) are V + (x) = — 2msn 2 (y + 
K(m)), — 2msn 2 (y + iK'(m)), — 2msn 2 (y + K(m) + iK'[m)) which are all strictly isospectral potentials to 
the original Lame potential. Thus, in this case we do not obtain any new solvable potentials by using 
supersymmetry. 

Now consider the case a = 2. All the five band edge eigenvalues and eigenfunctions of the PT- invariant 
Lame potential V{x) = — 6msn 2 (y) have already been given by us in Table 4 of ref. Starting from 

any of the five band edge eigenfunctions and calculating the corresponding superpotentials, we obtain five 
different supersymmetric partner potentials all of which have the same band edge energy eigenvalues as 
given in Table 4 of ref. In Table 5 we have given the expressions for these five different strictly 

isospectral potentials. It is worth noting that out of these five potentials, three are self-isospectral - 
they are the PT-invariant GAL potentials [2,2,2,0]. Hence, truly speaking, we only have three genuinely 
different potentials, all having the same band edge energies. For each of these cases, using the formalism 
of supersymmetric quantum mechanics [S], we can easily obtain expressions for the corresponding five 
eigenstates. Now, again by starting from these eigenfunctions, we can construct still different partner 
potentials but with identical band edges. In this way, one could construct a large number of periodic 
potentials with five band edges and two band gaps, all strictly isospectral to the PT-invariant Lame 
potential (jUJ with a = 2. 

Similarly, if we consider the PT invariant Lame potential ()44|) with a = 3, then we have 7 band edge 
eigenfunctions and eigenvalues all of which are analytically known and are given in Table 1 of ref. [I]. 
Again, using supersymmetry, we can obtain seven different partner potentials V+ all with the same band 
edge eigenvalues. By starting from any one of them and using other eigenfunctions recursively, we can in 
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principle construct a huge class of new isospectral potentials. Particular mention may be made of the case 
when we start from the eigenfunction sn(y)cn(y)dn(y) of the potential V(x) = — 12msn 2 (y). It is easily 
shown that the corresponding partner potential V + (up to a constant) is given by 

V(x) = -m[6sn 2 (y) + 2sn 2 (y + K(m)) + 2sn 2 (y + iK'(m)) + 2sn 2 (y + K(m) + iK'{m))\ . (45) 

Thus, we see that the PT-invariant GAL potential [6, 2, 2, 2] has precisely three bands, three band gaps 
and seven band edges, since it is the supersymmetric partner of the PT-invariant Lame potential (|44j) with 
a = 3. The process described above is readily extended to any Lame potential with integer a. We can 
start from any of the 2a + 1 band edges and obtain the corresponding supersymmetric partner potentials 
all having the same band edges. 

We have shown that the SUSY partners of the PT-invariant Lame potentials [6, 0, 0, 0] and [12, 0, 0, 0] 
are the potentials [2, 2, 2, 0] and [6, 2, 2, 2] respectively. What about the higher Lame potentials? In this 
connection, it is amusing to notice that the band edges of the PT-invariant Lame potential [20, 0, 0, 0] and 
the potential [6, 6, 6, 2] (which follow from Table 4) are identical. For example, out of 9 band edges, the 6 
band edge energy eigenvalues of [20, 0, 0, 0] are given by 

E = -5(m+2)±V4m 2 - 9m + 9 , E = -5(l+m)±2 v / 4m 2 + m + 4 , E = 5(l + 2m)±2\/9m 2 - 9m + 4 . 

(46) 

It is easily seen from Table 4 that exactly the same eigenvalues are obtained when a, 6, /, g take the 
values (2, 2, —3, 1), (2, —3, 2, 1), (—3, 2, 2, 1). Similarly, one can show that the three remaining eigenvalues 
of [20, 0, 0, 0] satisfy the same cubic equation as [6, 6, 6, 2] when a, b, f, g take the values (2, 2, 2, —2). 

In fact, one can show that the number (and structure) of band edges of the PT-invariant Lame potential 
[2a(2a+l), 0, 0, 0] is same as the QES states of the potential [a(a+l),a(a+l),a(a+l),(a—T)a\. For example, 
for this PT-invariant Lame potential it is well known that out of the 4a+l band edges of the Lame potential, 
a states each are of the form cn(y)sn(y)F a _i(sn 2 (y)), cn(y)dn(y)F a _i(sn 2 (y)), dn(y)sn(y)F a _i(sn 2 (y)), 
while the remaining a + 1 states are of the form F a (sn 2 (y)). Using Table 4, it is easily shown that there are 
again 4a + 1 QES states of the potential [a(a + 1), a(a + 1), a(a + 1), (a — l)a], out of which a QES states 
each are obtained when a, b, /, g are of the form a, a, a — 1, a — 1, or a, a — 1, a, a — 1, or a — 1, a, a, a — 1, 
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while a + 1 QES states are obtained when a, b, /, g are of the form a, a, a, —a. In fact we believe that all 
the band edge eigenvalues of the potentials [2a(2a + 1), 0, 0, 0] and [a(a + 1), a(a + 1), a(a + 1), (a — l)a] are 
identical. While this is easily shown for low values of a, at the moment, a general proof is still lacking. 

Similarly, one can show that the number (as well as the structure) of band edges of the PT-invariant 
Lame potential [(2a — l)2a, 0, 0, 0] is the same as the QES states of the potential [a(a + 1), (a — l)a, (a — 
l)a, (a — l)a]. For example, it is well known that out of the 4a — 1 band edges of the PT-invariant Lame 
potential, a states each are of the form cn(y)F a _i(sn 2 (y)), dn(y)F a _i(sn 2 (y)), sn(y)F a _i(sn 2 (y)), while the 
remaining a — 1 states are of the form sn(y)cn(y)dn(y)i ? a _2(sn 2 (y). Using Table 4, it is easily shown that 
there are 4a — 1 QES states of the potential [a(a + 1), (a — l)a, (a — l)a, (a — l)a], out of which a QES states 
each are obtained when a, b, f, g are of the form a, —a, a — 1, a — 1, or a, a — 1, —a, a — 1, or a, a — 1, a — 1, —a, 
while a — 1 QES states are obtained when a, b, f, g are of the form —a — 1, a — 1, a — 1, a — 1. In fact we believe 
that all the band edge eigenvalues of the potentials [(2a— l)2a, 0, 0, 0] and [a(a+l), (a— l)a, (a— l)a, (a— l)a] 
are identical. While this is easily shown for low values of a, a general proof is not available. 

On the basis of these results, we then conjecture that the potentials [a(a + l), a(a + l), a(a + l), (a — l)a], 
for integer a, have the same band edges as the Lame potential [2a(2a + 1), 0, 0, 0] and hence these potentials 
also have precisely 2a band gaps and (4a + 1) band edges, all of which are known in principle. Further, 
the potentials [a(a + 1), (a — l)a, (a — l)a, (a — l)a] have the same band edges as the Lame potential 
[(2a — l)2a, 0, 0, 0] and hence are also potentials with a finite number (2a — 1) of band gaps. It would be 
nice to have a general proof. 

3.2 Supersymmetry Partners of PT-invariant Associated Lame Potentials 

We start our discussion with the a = 2, b = 1 associated Lame potential and its corresponding PT-invariant 
potential V PT (x) = — 6msn 2 (y) — 2mcn 2 (y)/dn 2 (y) . All five band-edge eigenvalues and eigenfunctions for 
this potential have been given by us in Table 3 of ref. [IT] . As established previously [S1E!> this is a 
self-isospectral potential and hence using the band edge eigenfunction dn 2 (y) does not give any new partner 
potential. However, if instead we use the remaining four band edge eigenfunctions, then one gets four new 
SUSY partner potentials which are strictly isospectral to the PT-invariant [6,2,0,0] potential. 
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Let us now consider the PT-invariant AL potential [a(a+l), (a— 2)(a— 1), 0, 0], i.e. the potential (JSJ) with 
b = a — 2, / = g = 0. As shown by us jH], one of its exact band edge eigenfunction is ip(x) = cn(y)dn a_1 (y). 
It is easy to see that the corresponding partner potential V+ (up to a constant) is the potential [(a — 
l)a, (a — l)a, 2, 0]. Thus we immediately conclude that the PT-invariant potential [(a — l)a, (a — l)a, 2, 0] is 
strictly isospectral to the PT-invariant AL potential [a(a + 1), (a — 2) (a — 1), 0, 0]. In the special case when 
both a, b are integers, in view of our results on AL potentials [S] , it then follows that the GAL potential 
[(a — l)a, (a — l)a, 2, 0] has a band gaps and a bands, out of which b = a — 2 bands are rather unusual. 

Note that if instead we use ip(x) = sn(y)dn a_1 (y, m), which is also one of the exact eigenfunctions of 
the above AL potential, then nothing new is obtained. In particular, the corresponding partner potential 
is [(a — l)a, (a — l)a, 0, 2] which is strictly isospectral to the potential [(a — l)a, (a — l)a, 2, 0]. 

Let us now consider the PT-invariant AL potential [a{a + 1), (a — 3)(a — 2), 0,0], i.e. the potential 
(j3J) with b = a — 3, f = g = 0. As shown by us [Sj, one of its exact band edge eigenfunction is ip(x) = 
sn(y)cn(y)dn a ~ 2 (y). It is easy to see that the corresponding partner potential V+ (up to a constant) 
is the potential [(a — l)a,(a — 2)(a — 1),2,2] in the notation of ©• Thus we immediately conclude 
that when a, b are integers, then this PT-invariant potential is strictly isospectral to the AL potential 
[a(a + 1) , (a — 3) (a — 2) , 0, 0] , has a band gaps and a bands, out of which b = a — 3 bands are rather unusual. 

We can generalize the above arguments. In particular, we find that the number (and even structure) 
of the potentials [(a — p)(a — p + 1), (a — p — l)(a — p),p(p + l),p(p + 1)] is the same as the AL potentials 
[a(a + 1), (a — 2p — l)(a — 2p), 0, 0]. For example, as remarked in the previous section, if b = a — 2p — l(p = 
0, 1, 2, ...), then there are p eigenstates of the form sn(y)cn(y)dn(y)F p _i(sn 2 (y)), p + 1 eigenstates of the 
form dn a_2p (y)F p (sn 2 (y)), a—p eigenstates of the form cn(y)dn 2p+1-a (y)F a _ p _i(sn 2 (y)) and also a — p 
eigenstates of the form sn(y)dn 2p+1_a (y)i ? a _p_i(sn 2 (y)). Using Table 4 it is easy to show that for the GAL 
potential, [(a — p)(a — p + 1), (a — p — l)(a — p),p(p + l),p(p + 1)], there are p eigenstates of the form 
sn~ p (y)cn~ p (y) dn 1+p-a (y)F p _i(sn 2 (y)), p + 1 eigenstates of the form dn a_p (y)cn -p (y) sn~ p (y)F p (sn 2 (y)), 
a—p eigenstates of the form cn p+1 (y)sn -p (y) dn p+1_a (y) i ? a _ p „i(sn 2 (?/)) and also a — p eigenstates of the 
form sn p+1 (y)cn -p (y) dn p+1_a (y)i ? a _ p _i(sn 2 (y)). In fact we believe that all the band edge eigenvalues of 
the potentials [a(a + 1), (a — 2p)(a — 2p+ 1),0, 0] and [(a — p)(a — p + 1), (a — p — l)(a — p) , p(p + 1) , p(p + 1)] 
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are identical. While this is easily shown for low values of a and p, a general proof is still lacking. 

Similarly, we can show that the number (and even structure) of the potentials [(a — p)(a — p + 1), (a — 
p)(a — p + l),p(p + l),(p — l)p] is the same as the AL potentials [a(a + 1), (a — 2p)(a — 2p + 1),0,0]. 
In particular, for the AL potential, as shown in Sec. 2, when b = a — 2p, then there are p eigen- 
states of the form cn(x)dn a_2p+1 (x)-F p _i(sn 2 (x)), p eigenstates of the form sn(x)dn a_2p+1 (x)i ? p _i(sn 2 (x)), 
a — p eigenstates of the form sn(x)cn(x)dn 2p ~ a (x)F a _ p _i(sn 2 (a;)), and a — p + 1 eigenstates of the form 
dn 2p_a (x)F a _ p (sn 2 (x)). It is easily shown that for the potential [(a — p)(a — p+ 1), (a — p)(a — p + l),p(p + 
1), (p — l)p], there are 4a — 1 QES states of similar form. In particular, there are p eigenstates of the form 
sn~ p (y)cn 1_p (y) dn p ~ a (y)F p _i(sn 2 (y)), p eigenstates of the form dn a+1 ~ p (y)cn 1_p (y) sn~ p (y)F p _i(sn 2 (y)), 
a — p + 1 eigenstates of the form cn p (y)sn~ p (y) dn p_a (y)F a _ p (sn 2 (y)) and also a — p eigenstates of the form 
sn p+1 (y)cn 1_p (y) dn p ~ a (y) i ? a _ p _i(sn 2 (y)). In fact we believe that all the band edge eigenvalues of the 
potentials [a(a + 1), (a — 2p)(a — 2p + 1), 0, 0] and [(a — p)(a — p + 1), (a — p)(a — p + l),p(p + 1), (p — l)p] 
are identical. While this is easily shown for low values of a and p, we don't yet have a general proof. 

3.3 SUSY Partners of Potentials with b = f = 

Let us now consider the SUSY partners of the potential [a(a + 1), 0, 0, g(g + 1)] which for integral values of 
a, g, is a problem with a finite number of band gaps. By exactly following the above discussion about the 
PT-invariant AL potential, we can construct a host of new potentials with a finite number of band gaps. 
For example, by starting from the potential [6, 0, 0, 2] and following the procedure as in the AL case, we 
can easily obtain four new SUSY partner potentials, all with two band gaps. 

From Table 1 we observe that for integral a, two of the exact eigenfunctions of the potential [a(a + 
1),0, 0, (a — 2)(a — 1)] with a band gaps are cn(y)sn a ~ 2 (y) and dn(y)sn a ~ 2 (y). It is easily seen that if 
we start with either of these eigenfunctions, then the corresponding SUSY partner potential with the 
same finite (a) number of band gaps is the potential [(a — l)a, 2, , 0, (a — l)a] (or its isospectral partner 
[(a-l)a,0,2,(o-l)a]). 

From Table 1 we also observe that one of the exact eigenfunction of the potential [a(a + 1),0,0, (a — 
3) (a — 2)] is cn(y, m)dn(y, m)sn a ~ 3 (y, m). On starting with this eigenfunction, it is easily shown that the 
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corresponding SUSY partner potential is [(a — l)ct, 2, 2, (a — 2)(a — 1)] which therefore must also be a 
potential with finite (a) number of band-gaps in case a is an integer. 

Similarly, by starting from the finite band-gap potentials [a(a + 1), 0, 0, (a — 2p — l)(a — 2p)\ as well as 
[a(a + 1), 0, 0, (a — 2p)(a — 2p + 1)], and following the discussion in the case of PT-invariant AL potential, 
it is easily shown that the corresponding SUSY partners with the same (finite) number of band gaps are 
the potentials [(a — p)(a — p + l),p(p + l),p(p + l),(a — p— l)(a — p)] and [(a — p)(a — p+l),p(p + l), (p - 
l)p, (a — p)(a — p + 1)] and respectively, where a and p are positive integers. 

3.4 SUSY Partners of Potentials with b = g = 

Let us now consider the SUSY partners of the PT-invariant potential [a(a + 1),0, /(/ + 1),0]. From 
Table 2 we observe that two of the exact eigenfunctions of the potential [a(a + 1), 0, (a — 2) (a — 1), 0] arc 
sn(y)cn a ~ 2 (y) and dn(y)cn a_2 (y). It is easily seen that if we start with either of these eigenfunctions, then 
the corresponding SUSY partner potentials turn out to be [(a — l)a, 2, (a — l)a, 0] or [(a — l)a, 0, (a — l)a, 2]. 
Since we know that the potentials [a(a + 1), a(a — 1), 2, 0] as well as [a(a + 1), 2, 0, (a — l)a] have a finite 
number of band gaps, we conjecture that maybe the potential [a(a + 1), 2, (a — l)a, 0] also has only a finite 
number (a) of band gaps when a is an integer. 

From Table 2 we also observe that one of the exact eigenfunctions of the potential [a(a + 1),0, (a — 
3) (a— 2), 0] is sn(y)dn(y)cn a ~ 3 (y). Starting with this eigenfunction, it is easily shown that the corresponding 
SUSY partner potential is [(a — l)a, 2, (a — 2)(a — 1),2]. Again, since for integer a, the potential [(a — 
l)a, (a — 2) (a — 1), 2, 2] has only a finite number of band gaps, it is tempting to conjecture that the same 
may also be true for the potential [(a — l)a, 2, (a — 2) (a — 1), 2]. 

Similarly, by starting from the finite band-gap potentials [a(a + 1), 0, (a — 2p — l)(a — 2p), 0] as well as 
[a(a + 1), 0, (a — 2p)(a — 2p+ 1), 0], and following the discussion in the case of PT-invariant AL potentials, 
it is easily shown that the corresponding SUSY partners with the same number of band gaps are the GAL 
potentials [(a — p)(a — p + l),p(p + + 1), (a — p — l)(a — p)\ and [(a — p)(a — p + + 1), (p — 

(a — p)(a — p + 1)] respectively when a and p are integers. 
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3.5 SUSY Partners of Potentials with / = 

Let us now consider the SUSY partners of the potential [a(a + 1) , 6(6 + 1), 0,g(g + 1)]. From Table 3 we 
observe that one of the exact eigenfunctions is dn _fe (y)sn~ 9 (y) when a + 6 + g = 0. If we start with this 
eigenfunction, then the corresponding SUSY partner potential turns out to be [(a — l)a, (6 — 1)6, 0, (g — l)g]. 

From Table 3 we also observe that an exact eigenfunction of the potential [a(a + 1), 6(6+ 1), 0,g(g + 1)] 
is cn(y)dn _6 (y)sn -9 (y) when a + b + g = 1. Starting with this eigenfunction, it is easily shown that the 
corresponding SUSY partner potential is [(a — l)a, (6 — 1)6, 2, (g — l)g]. 

In summary, we have discovered a large number of complex PT-invariant periodic potentials with a 
finite number of band gaps, many occurring when the parameters a, 6, c, d have specific integer values. This 
leads us to make the plausible conjecture that all GAL potentials (j3J) for integer values of a, 6, /, g have a 
finite number of band-gaps, but there is as yet no formal proof. 

4 Heun's Equation and the Generalized Associated Lame Equation 

In this section, we point out an interesting connection between Heun's differential equation [7J and the 
generalized associated Lame equation (JSJ) . This connection enables us to use the various solutions of eq. 
(jSJ) obtained in this paper to write down several solutions of Heun's equation which have apparently not 
been studied in the mathematics literature. 

The canonical form of Heun's equation is given by [7] 

d 2 7 5 e d af3x — q 



G(x) = 0, (47) 



_dx 2 x x — 1 x — c dx x(x — l)(x — c) 

where a, (3, 7, 5, e, q, c are real parameters, except that c 7^ 0, 1 and the first five parameters are constrained 
by the relation 

7 + «5 + e = Q + /3 + l. (48) 
If we make the transformation x = sn 2 (y,m), then Heun's equation takes the form [7J 

F "{y) + [(i - 2e)m ^M + (1 _ 26 f_^y}Myl + (27 _ if-^yl ]F ' {y) 

dn(y) cn(y) sn(y) 

-[Amq - 4a(3msn 2 (y)]F(y) = , (49) 
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where [G(x) = F(y)] and m = 1/c. It is interesting to note that eq. (|49[) is very similar to the 4> equation 
(|14jl which we have analyzed in great detail. In particular, with the identification 

b=~-e, / = i - J , g = ~ - 7 , 6 + / + 5 = ^- a-/?,4a/3 = Q,4mg = J R, (50) 

all the results discussed above can be immediately used to obtain different solutions of Heun's equation. 
It turns out that using the mid-band states obtained in Sec. 2, one generates new quasi-periodic solutions 
of Heun's eq. (|49|) . which we discuss in a separate publication [S]. 
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Table 1: Energy eigenstates of PT-invariant GAL potentials with b = f = 0, g = n — a, n = 0, 1, 2, . . .; 

5i = + m) 2 (a - l) 2 - (2a - l)(2a - 3)m , <5 2 = V[a - 1 + m(a - 2)] 2 - (2a - l)(2a - 5)m , 
<5 3 = v/[a - 2 + m(a - 1)J 2 - (2a - l)(2a - 5)m , <5 4 = ^(1 + m) 2 (a - 2) 2 - (2a - l)(2a - 7)m . 



n 5(5 + 1) £ sn- a (y)V> 








(a — l)a 






-(1 + m)a 2 


1 
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1 


(a 


-2) (a - 


1) 




-a 2 - m(a - l) 2 


cn(j/) 
sn(j/) 


1 


(a 


-2)(a- 


1) 




-(a - l) 2 - ma 2 


dn(j/) 
sn(j/) 


2 


(a 


-3)(a- 


2) 




-(l + m)(a- l) 2 


cn(j/)dn(j/) 
sn 2 (y) 
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(a 


-3)(a- 


2) 




-(1 + m)(a 2 - 2a + 2) ± 2<5i 


[£+(l+m)(a-2) 2 ]sn 2 (j/)+2(2a-3) 
sn 2 (j/) 
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(a 


-4)(a- 


3) 


"(a 2 


- 2a + 2) - (a 2 - 4a + 5)m ± 25 2 


[[£+(a-2) 2 +m(a-3) 2 ]sn 2 (j/)+2(2a-5)]cn(j/) 
sn 3 (y) 


3 


(a 


-4)(a- 


3) 


"(a 2 


- 4a + 5) - (a 2 - 2a + 2)m ± 25 3 


[[£+(a-3) 2 +m(a-2) 2 ]sn 2 (j/)+2(2a-5)]dn(j/) 


sn 3 (y) 


4 


(a 


-5)(a- 


4) 




-(l + m)(a 2 -4a + 5)±2<5 4 


[[£+(l+m)(a-3) 2 ]sn 2 (j/)+2(2a-7)]cn( J/ )dn( J/ ) 




sn 4 (j/) 



Table 2: Energy eigenstates of PT-invariant GAL potentials with parameters 6 = 5 = 0, / = n — a, 
n = 0,1,2,...; 05 = v/(a- 1 + m) 2 - (2a - l)m , 5 6 = y/(a - 1 + 2m) 2 - 3(2a - l)m , 
<5 7 = v 7 ^ - 2 + 2m) 2 - (2a - l)m , <5 8 = y/(a - 2 + 3m) 2 - 3(2a - l)m . 
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Table 3: Energy eigenstates of PT-invariant GAL potentials with / = 0, g = n — a — b, n = 0, 1, 2, . . .; 

S 9 = + m)(a - 1) + b} 2 - (2a - l)(2a + 2b- 3)m, 
<5io = v/[a + b - 1 + m(a - 2)J 2 - (2a - l)(2o + 26 - 5)m. 



n 



E dn 6 (y)sn-( a+6 )(y)^ 

-(a + 6) 2 -ma 2 1 

1 _ (a + 6) 2_ m(a _ 1)2 2g 

2 -(l+m)-(a + 6-l) 2 -m(a-l) 2 ±2«y 9 ^ + (a + b- 2 ) 2+m (aJ)W(,) + 2(2a + 2b-3) 

3 _(l + m )_ (a + 6 _l)2_ m(a _ 2) 2 ±2(yi0 [[g+(a + 6-2)'+m( a -3^(y)+2(2a+2 & -5)lcn( y ) 



Table 4: Energy eigenstates of PT-invariant GAL potentials with / = 2n — a — b — g, n = 0, 1, 2, . . .; 
tfii = y/[{a + b-l)+m(l-b- g)] 2 - (2a - 1)(1 - 2g)m . 



n 


E 


sn3(y)dn b (y)cn-( a+b +9) (y)^ 



1 
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Table 5: The five supersymmetric partner potentials of the PT-invariant Lame potential VE T (x) = 
— 6rasn 2 (y) . Here y = ix + (3 and 5 = y/l — m + m 2 . All partner potentials have a period 2K'(m). 
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